Universal and non-universal tails of distribution functions 
in the directed polymer and Kardar-Parisi-Zhang problems 
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The optimal-fluctuation approach is applied to study the most distant (non-universal) tails of 
the free-energy distribution function Pl{F) for an elastic string (of a large but flnite length L) 
interacting with a quenched random potential. A further modification of this approach is proposed 
which takes into account the renormalization effects and allows one to study the closest (universal) 
parts of the tails. The problem is analyzed for different dimensions of a space in which the polymer 
is imbedded. In terms of the stochastic growth problem, the same distribution function describes the 
distribution of heights in the regime of a non-stationary growth in the situation when an interface 
starts to grow from a flat conflguration. 
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I. INTRODUCTION 

A large variety of physical systems can be described 
in terms of an elastic string interacting with a quenched 
random potential. The role of such a string can be played 
by a domain wall in a two-dimensional magnet, a vortex 
line in a superconductor, a dislocation in a crystal, and 
so on; however following Ref . HI the systems of such a kind 
are usually discussed under the generic name of a directed 
polymer in a random medium. The unfading interest to 
this problem is additionally supported by its resemblance 
to more complex systems with quenched disorder (e.g., 
spin glasses), as well as by its close relation to the dy- 
namics of a randomly stirred fluid and to the problem of 
a stochastic growth (see Refs. andlH for reviews). 

One of the main objects of interest in the directed poly- 
mer problem is Pl (F) , the free-energy distribution func- 
tion for large polymer length L. In particular, the knowl- 
edge of this distribution function allows one to make con- 
clusions on the distribution of displacements. The first 
important step in the analysis of Pl (F) was made twenty 
years ago by Kardar^i who suggested that all moments of 
Pl (F) can be found by calculating the moments Z„ = Z"- 
of the distribution of the partition function Z and pro- 
posed an asymptotically exact method for the calculation 
of Zn in a (1 + 1) -dimensional system (a string confined 
to a plane) with a (5-correlated random potential. How- 
ever, soon after that Medina and Kardar^ understood 
(see also Ref. that the information provided by the 
approach introduced in Ref. 4 is insufficient for finding 
any of the moments of Pl{F). However, it allows one to 
find^ the tail of Pl (F) at large negative F (the left tail) . 
In such a situation the conclusions on the width of the 
distribution function have to rely on the assumption that 
at large L it acquires a universal form. 
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incorporating the dependence on all parameters through 
a single characteristic free-energy scale F^,{L) cx L'^, 



which therefore can be extracted from the known form 
of the tail. The form of Eq. ([1]) assumes that F, the free 
energy of a directed polymer in a given realization of a 
disorder, is counted off from its average, F{L) or, more 
precisely, from the linear in L contribution to F (that 
is, L\imL-*oo[F{L)/ L]). The same convention is implied 
below. 

Only recently it has been understood^ that the form of 
the tail following from Zhang's analysis^ is applicable for 
the description only of the most distant part of the left 
tail and therefore has no direct relation to the universal 
form of the distribution function which is achieved in the 
limit oi L ^ OO. At large but finite L the form of Pl{F) 
given by Eq. ([1]) can be expected to be achieved only for 
not too large fluctuations of F [that is, for |F| <C Fc{L) 
with Fc{L)/Fif{L) tending to infinity with the increase 
of L], whereas the behavior of Pl{F) at |F| 3> Fc{L) re- 
mains nonuniversal and is not obliged to have anything in 
common with P^{F / F^,) . In particular, it can incorporate 
quite different characteristic free-energy scales. Thus, the 
fact that Zhang's approacb^^ reproduces both the correct 
form of the left tail of (F/ F^, ) and the correct estimate 
of the universal free energy scale {L) (which is the only 
relevant free-energy scale inside the universal region), is 
not more than a happy coincidence. In contrast to that, 
the behavior of the right tail of Pl (F) inside the univer- 
sal region is qualitatively different from its behavior in 
the nonuniversal part of the tailf^ 

In this article the analysis of the universal and non- 
universal tails of Pl (F) developed in Ref. is presented 
in more detail and also is extended to the investigation 
of (1 -|- c?) -dimensional systems, in which polymer's dis- 
placement can be treated as a c?-dimensional vector. The 
article is organized as follows. 

In Sec. |TT] we formulate the continuous model which 
is traditionally applied for the description of the di- 
rected polymer problem and briefly review its relation 
to the Kardar-Parisi-Zhang (KPZ) modelS of a stochas- 
tic growth, as well as to the Burgers turbulence problem. 

Sec. [HT] provides a short introduction to the optimal 
fluctuation approach, which can be used for the descrip- 
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tion of the most distant (non-universal) parts of the tails 
of Pl{F). In Refs. [l^ an analogous approach has been 
used to investigate the distribution of velocity and its 
derivatives in the Burgers turbulence problem, which 
however requires one to consider optimal fluctuations 
with completely different structures then studied here. 

In Sec. IIVI the optimal fluctuation approach is applied 
for the analysis of the far- left tail of P]^{F), and in SecfVl 
of the far-right tail. Our main attention is focused on the 
systems with a (5-correlated random potential; however 
for d > 2 the problem with purely J-functional corre- 
lations becomes ill-defined, so we also consider the case 
when a random potential correlations can be character- 
ized by a finite correlation radius. 

For finding the universal parts of both tails one also 
has to look for optimal fluctuations, but taking into ac- 
count that in this regime the parameters of the system 
have to be considered as scale dependent due to their 
renormalization by fluctuations. This is done in Sec. IVII 
The validity of this approach is confirmed by the con- 
sistency of its predictions with the results of the exact 
solutionii of the (1 + 1) -dimensional polynuclear growth 
(PNG) model, as well as by obtaining identical estimates 
for -F»(-L) in the left and right tails. 

The concluding Sec. lVIIl is devoted to summarizing the 
results and comparing them with some results of other 
authors, whereas in Appendix 1X1 we discuss how some of 
the results of this work can be derived in terms of the 
Kardar- Zhang replica approach 4*^ 

Our main attention throughout this work is focused on 
a system with free initial condition, that is, we assume 
that only one end of a string is fixed, whereas the other 
one is free to fluctuate. In terms of the KPZ problem^ 
the same distribution function describes the distribution 
of heights in the regime of a nonstationary growth in 
the situation when an interface starts to grow from a flat 
configuration {L being the total time of the growth). One 
only has to bear in mind that the height (as defined in the 
standard form of the KPZ equation) and the free energy 
of the directed polymer problem differ from each other 
by the sign. Therefore, what we call here the left (right) 
tail of Pl{F) in terms of the KPZ problem corresponds 
to the right (left) tail of the height distribution function. 

Finally, Appendix |B] is devoted to demonstrating that 
when both end points of a directed polymer are fixed, the 
form of the left tail of Pl (F) remains basically the same 
as for free initial condition. 



II. THE MODEL 

We consider an elastic string in a (1 -I- d) -dimensional 
space interacting with a random potential V{t,x). The 
coordinate along the average direction of the string is 
denoted t for the reasons which will become evident few 
lines below. Such a string can be described by the Hamil- 



tonian, 
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where the first term describes the elastic energy and the 
second one the interaction with a random potential. Note 
that the form of the first term in Eq. ^ relies on the 
smallness of the angle between the string and its preferred 
direction. 

The partition function of a string which starts ai t — 
and ends at the point {t, x) is then given by the functional 
integral. 
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x(t)=x 



x(0)=x' 



z(t,x)= / dx'z(0,x') / X'x(t') exp (-i//T) , 

(3) 

where T is the temperature. Naturally, z{t,x.) depends 
on the initial condition ai t — 0. The fixed initial condi- 
tion, x(< = 0) = Xq, corresponds to z{0,x) — S{x — xq), 
whereas the free initial condition (which implies the ab- 
sence of any restrictions on x at t = 0) to 



z{0, x) = const . 



(4) 



Since Eq. ([3]) has exactly the same form as the Eu- 
clidean functional integral describing the motion of a 
quantum particle whose mass is given by J in a time- 
dependent random potential V{t,x) (with t playing the 
role of imaginary time and T - of Plank's constant h), 
the evolution of z(t,x) with the increase in t has to be 
governed by the imaginary-time Schrodinger equation 



T 



dz 



2J 



V{t,x) 



z{t,x) . 



(5) 



As a consequence of this, the evolution of the free energy 
corresponding to z{t,x), 



/(t,x) = -Tln[z(i,x)] 



(6) 



is governedi^ by the Kardar-Parisi-Zhang (KPZ) 
equation;^ 



1 



df 



dt 2J 



(7) 



with the inverted sign of /, where t plays the role of time 
and 1/ = T/2J of viscosity. On the other hand, the deriva- 
tion of Eq. ([7]) with respect to x allows one to establish 
the equivalence^^ between the directed polymer problem 
and the Burgers equationi^ with random potential force. 



du 1 

at ^ 2 

where the vector 



Vu^-j/V^u= -VV"(t,x), (8) 



U(t,x) EE jV/(i,x) 



(9) 
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plays the role of velocity. Note that in terms of the 
KPZ problem the free initial condition ^ corresponds to 
starting the growth from a flat interface, /(0,x) = const, 
and in terms of the Burgers problem to starting the evo- 
lution from a liquid at rest, u(0,x) — 0. 

To simplify an analytical treatment, the statistic of a 
random potential is usually assumed to be Gaussian with 



V{t,x)^0, V{t,x)V{t',yi') = S{t~t')U{x~x') , (10) 

where an overbar denotes the average with respect to dis- 
order. Our main attention below is focused on the case 
of purely (5-functional correlations, U{x.) — UoS{x.). How- 
ever, for d > 2 the problem with such a form of cor- 
relations is ill-defined and needs a regularization, so we 
also consider the case when C/(x) can be characterized 
by a finite correlation radius ^. On the other hand, we 
always assume that the correlations in the t direction are 
(5-functional, because in almost all situations considered 
below the finiteness of the correlation radius in the t di- 
rection can be ignored as soon as it is small in comparison 
with the total length of a string. 



III. OPTIMAL-FLUCTUATION APPROACH 

When the distribution of V{t, x) is Gaussian and satis- 
fies Eqs. PH)) . the probability of any realization of V{t, x) 
is proportional to exp[— where the action S{V} is 
given by the functional 



S{V} 



dt 



dxdx'V(t,x.)U^\x-x')V{t,yi') 



(11) 

Here U ^(x) denotes the function whose convolution 
with [/(x) is equal to i5(x). Accordingly, the probabil- 
ity of any time evolution of /(t,x) is determined by the 
action S{f}, which is obtained by replacing V(t,x) in 
Eq. ini) by the left-hand side of the KPZ equation 0. 

To find the most optimal fluctuation having the largest 
probability (in literature it is often called "instanton"), 
one has to minimize S{f} for the given boundary condi- 
tions at t — and t = L. A convenient way to perform 
such a minimization consists in replacing S{f} by 



S{f.tA^j^ dt^j dy 



1 



df 



dt 2J 



(V/)2 - i.VV 



dxdx'/^(t,x)C/(x-x')A*(i,x')| (12) 

where fi{t,x) is an auxiliary field with respect to which 
S{f,fi} also has to be extremized. Variation of Eq. IT^ 
with respect to fi{t, x) reproduces the KPZ equation ([T]) 
with 



V{t,x)^ J dx'J7(x-x')M(i,x') 



(13) 



whereas its variation with respect to f{t,x.) leads to 

dfi/dt + div{un)+i^\/^fi = 0, (14) 



where u(t,x) = V f{t,x.)/ J is the "velocity" entering the 
Burgers equation ([8]). The form of Eq. (|14p implies that 
the integral of fi{t,x) over whole space is a conserved 
quantity, whereas substitution of Eq. into Eq. ([TT|) 
shows that in terms of fi(t, x) the action can be rewritten 
as 

S-j/ij^i J '^^ J J d:>idx' fi{t,x)U{x^ x')fi{t,x') . 

(15) 

In a system with (5-functional correlations, 
U{x) = Uq5{x), V and /z differ from each other 
only by a constant factor Uq, and accordingly, Eq. ^ 
can be replaced by 



dV/dt + div(uy) + i/V V = . 



(16) 



If the beginning of a polymer (at i = 0) is not fastened 
to a particular point and is free to fluctuate, the initial 
condition for the partition function z{t, x) has to be cho- 
sen in the form z{0,x) = const. In such a case to find 
the tails of Pl{F) one has to find the solution of Eqs. ([7]) 
and (11411 which satisfies the initial condition 



(17) 



/(0,x) = 0, 
and the final condition 



(18) 



where for the left tail F < and for the right tail F > 
0. Alternatively, condition can be imposed by the 
inclusion of the (5-functional factor, 



/ 



dAexp{fA[/(L,x = 0)-F]}, 



(19) 



into the functional integral defining the probability of 
a fluctuation. In such a case condition (fT8|) for /(i,x) 
should be replaced by the condition for /i(i,x). 



^l{L,x) = HoS{x) , 



(20) 



where, however, the value of /io oc A has to be chosen to 
satisfy Eq. (fTS]) . 



IV. FAR-LEFT TAIL 

It turns out that in the case of the left tail the solution 
of Eqs. ^ and (fT4|) which satisfies boundary conditions 
(fT7|) and (|18|) can be constructed on the basis of the so- 
lution of these equations in which the potential V and all 
derivatives of / do not depend on t, which means that the 
time dependence of f{t,x) is decoupled from its spacial 
dependence and is as trivial as possible. 



f{t,x)^E{t-h) + fix), 



(21) 



where ti = const and E — const < 0. Below we for 
brevity call such solutions stationary. 
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For /(<,x) of form ((2T|) the replacement 
/(x) --rin*(x) 



(22) 



transforms the KPZ equation ([7]) into a stationary 
Schrodinger equation: 



(23) 



for a single-particle quantum- mechanical problem defined 
by the Hamiltonian 



H 



'2J 



nx), 



(24) 



where J plays the role of mass and T of Plank's con- 
stant h [compare with ([5])]. On the other hand, when 
both u = -(r/J)V^'/* and /x do not depend on t, Eq. 
(fT4|l is automatically fulfilled as soon as 



/i(x)(xvl/2(x), (25) 

which implies 

y(x) = -^y" f^x' f/(x - x')4'^(x') , (26) 

where A is an arbitrary constant. Substitution of Eq. ((26|) 
into Eq. ([^5]) allows one to replace them by a single non- 
linear Schrodinger equation, 



2 J 



/+(X} 
dx'C/(x-x')*^(x'). 
-oo 



Equation (|26p has been derived by Halperin and Lax^"* 
when looking for the optimal fluctuation of the potential 
y(x), which for the given value of the ground state en- 
ergy _B < of the quantum-mechanical Hamiltonian 
minimizes the functional 



s{V} 



dxdx'y(x)t/-i(x-x')V'(x'), (27) 



determining the probability of V{'x.) (or, equivalently, 
minimizes E for the given value of s{l^}). Apparently, in 
terms of our problem s{V} is related to the action 5'{V^} 
defined by Eq. pT|) as 5 = Ls. In the case of 5-functional 
correlations, C/(x) = C/o(5(x), and independent potential 
V{x), functional (fTTj) is reduced to 



A. 5-functional correlations, d = 1 



(28) 



In a 1 -|- 1-dimensional system with a (5-correlated ran- 
dom potential, U{x) = UoS{x), the localized solution of 
Eqs. (pS)) and (pS)) (the soliton) exists for any E < and 
can be found exactly jAi 



V{x) = 



-2E 
2E 



1/2 



cosh(a;/A) 



coah^ix/A) ' 



(29) 
(30) 



a) u(x) 



X 



b) 



f(x)i 




FIG. 1: The spacial dependence of u and / in the stationary 
solution of Eqs. (O and ^U^. 



where the length-scale 



(-2j£;)i/2 



(31) 



can be called soliton width. This allows one to conclude 
that the stationary solution of Eqs. ([7]) and (fT6|) is given 
by Eq. §U\) and 



f{t,x) = E{t - h) +T\n (2 cosh ^ 



(32) 



which follows from the substitution of Eq. into Eqs. 
([2T|) and ([22]). Note that in Eq. ^ the constant h has 
been redefined in order to absorb A. 

Differentiation of Eq. (|32|) with respect to x gives a 
stationary profile of u{x), 



u{x) — V tanh 



A 



schematically shown in Fig.[T]^a). Here 
V = T/JA 



(33) 



(34) 



is the velocity of the outward fiow created by the forces 
acting inside the soliton. The profile ([55]) up to a sign 
coincides with the one in a stationary shock wave with 
the same amplitude v. The solitons of such a kind (both 
stationary and moving ) have been discussed in a number 
of works by Fogedbyj^ 

The stationary profile of / described by Eq. ([5^ is 
schematically shown in Fig. [TJb). With the increase in 
time it is moving downward as a whole with a constant 
velocity df/dt = E. Away from the soliton's core, that 
is at I a; I 3> A, the dependence described by Eq. ([5^ can 
be approximated as 



fit,x)^ E{t-h) + i~2JE)^^^\x\. 



(35) 
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a) u(x) I 



1 



X 



b) f(x)t 




FIG. 2: The spacial dependence of u and / in the solution of 
Eqs. (O and JTe]) corresponding to the left tail of Pl{F). The 
arrows show the directions of motion of the two shock waves. 



Since Eq. ([55)1 describes a solution of the noiseless KPZ 
equation, its form does not depend on the form (or am- 
plitude) of the random potential correlator U{x). 

The stationary solution minimizes the action for the 
given negative value of df /dt = E. Therefore, it allows 
one to find the optimal value of S in situations when it 
is not influenced by the initial condition. Substitution of 
Eq. ((301) into Eq. ^ then gives 



5(A) 



2 T^L 



3 C/o J2A3 



(36) 



Apparently, the condition /(i, 0) — /(0, 0) = is fulfilled 
when E = F/L, which corresponds to 



F : 



rp2 



2JA2 



and 



S{F) 



Ay/2 T{-Ff/^ 
3 UoJ^/^L^/^ 



(37) 



(38) 



However, the real optimal fluctuation also has to re- 
spect the initial condition and it is clear that the spacial 
dependence of / in Eq. (j32|) in no way resembles the ini- 
tial condition (jl7p . In terms of the quantum- mechanical 
problem with time-independent potential V^(x) it is clear 
that the applicability of the relation F « EL requires 
to have {SE)L 3> T, where SE is the energy gap sepa- 
rating the ground state of the Hamiltonian ([M]) from the 
first excited state. Since in potential ([30)1 there exists only 
one bound level with a negative energy,>i£ whereas excited 
states can have any non-negative energy, this condition 
is equivalent to — F ^ T. 



For constructing a non-stationary solution which elim- 
inates the inconsistency between the forms of the sta- 
tionary solution and of the initial condition (without in- 
creasing the action), one has to complement the soliton 
shown in Fig. [ija) by two traveling shock waves [as shown 
in Fig.[2l^a)], whose existence does not require any addi- 
tional pumping. Both these shock waves will be moving 
outwards with velocity v/2. Their presence will change 
the profile of f{t,x) to the one shown in Fig. [UJb), so 
that f{x) will be given by Eq. ([32| only in the interval 
where f{t,x) < 0, whereas outside of this region it will 
coincide with the initial condition (|17p (with a smooth 
crossover between the two solutions) . This means that if 
a potential localized in the vicinity of x = is switched 
on at t = 0, its influence on f{t, x) a.t t > extends only 
to a finite (but growing with t) region, which is perfectly 
logical. 

In such a situation the constant ti in Eq. ((32)) [or in Eq. 
P5|) ] will have the meaning of an effective time required 
for the formation of the non-stationary solution shown 
in Fig. [51 At the initial stage, that is, at i < ti, the 
spacial distribution of V{x) will substantially differ from 
the one given by Eq. ([30]) The value of ti can be estimated 
from the comparison of the soliton width A = 2z//u with 
the velocity v of the flow it creates, which gives ti ^ 
vjv^ . This allows one to expect that for t\ the main 
contribution to the action is coming from the region in t 
where Eq. ([30]) gives a sufficiently accurate description of 
the solution, and therefore the value of S(F) is given by 
Eq. p8)) . In terms of F the constraint t\ corresponds 
to the condition —F ^ T (which was already derived 
above in different terms). 

The same condition allows one to neglect the final stage 
of the optimal-fiuctuation evolution. At this stage the 
potential has to shrink from the form given by Eq. (|30p 
to a 5-function as suggested by Eq. ([20| . Simultaneously, 
the downward tip of jix) has to change its shape from 
rounded to more sharp. The decrease in /(a; = 0) related 
to this process can be expected to be comparable with 
the change in / induced by the rounding of the tip, which 
according to Eq. ([32]) is of the order of T, and therefore 
for — F ^ T can be ignored. 

Note that the same answer for S{F), Eq. ([55]) . can 
be also obtained in the framework of the Kardar-Zhang 
replica approach based on mapping a system to a set 
of interacting bosons and keeping only the ground state 
contribution to the partition function of these bosons (see 
Appendix |X] for more details). 

In Appendix [B] we demonstrate that the change of the 
initial condition from free to fixed does not change the 
form of the main contribution to S{F). The same conclu- 
sion is even more easily attained in terms of the replica 
approach (see Appendix E| . 

In the remaining part of this section, we analyze the 
systems with an arbitrary dimension and/or finite-range 
correlations assuming that the main features of the opti- 
mal fluctuation determining the far-left tail are the same. 
Namely, we expect that in a growing region around x = 0, 
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the solution is close to the stationary solution, whereas 
outside of this region it is close to the initial condition 
f{t,x) = 0, the crossover between the two regions be- 
ing described by a corresponding solution of the noise- 
less KPZ equation. In such a situation the action of the 
optimal fluctuation is determined by the form of the sta- 
tionary solution. 



B. Generalization to d 7^ 1 

If the dimension of the transverse space d is not equal 
to 1, the joint solution of Eqs. (|23|) and ([26|) . that is, 
the wave function ^'(x) which minimizes the sum of a 
positive kinetic energy. 



IC 



T2 /d''x|V*(x)|2 
2J /d''x|1'(x)|2 



and a negative potential energy, 

/d''xT/(x)|*(x)|2 



V 



/d'ix|*(x)|2 



(39) 



(40) 



for a given value of the functional S"!!^} defined by Eq. 
(I28p . cannot be found exactly. However, in a situation 
when this wave function and, therefore, the potential 
V{x.) (X — [/o^'^(x) are well localized at some length scale 
A, an estimate for A and a qualitative relation between 
S and F can be obtained without finding the exact form 
of 4'(x). 

When ^'(x) can be characterized by a single relevant 
length-scale A, one has 



/C(A) 



JA2 ' 



(41) 



whereas the absolute value of V ^ ^(0) at a given S can 
be estimated with the help of Eq. (^5]) . which gives 



5^ Aav, 



and therefore 



V(A) 



SUq 



1/2 



(42) 



(43) 



For < d < 4 the sum /C(A) + V(A) has a minimum with 
respect to A when /C(A) ^ — V(A) and therefore both K, 
and —V have to be of the same order as —E = —F/L. 

Substitution oi V ^ F/L into Eq. (|42|l allows one to 
rewrite this relation as 



^dp2 
UnL 



(44) 



On the other hand, an estimate for A in terms of F can 
be obtained from the relation K, ^ —E, which gives 



T 

A(F)~- 



JL 



1/2 



(45) 



After that to obtain an estimate for S{F) one needs only 
to substitute Eq. ([i5|) into Eq. which leads to 



S{F) 



T'^(-F) 



2-d/2 



(46) 



Naturally, for d = 1 Eq. (|46| is consistent with Eq. ([SS]) 
derived in the Sec. lIV Al on the basis of the exact solution 
of Eqs. and 

For d > 4 the sum of /C(A) and V(A) at a given S is 
not bounded from below and tends to —00 when A 0. 
Accordingly, for any F < it becomes possible to find 
a stationary fluctuation with an arbitrary low action, so 
the method of optimal fluctuation is no longer applicable. 
However, it turns out that the range of the applicability of 
Eq. (|46p is even more narrow than the interval < d < 4, 
where the action of stationary fluctuations has a well- 
defined positive minimum. 

The point is that L enters Eq. (pS)) as the total time 
of the development of the optimal fluctuation of f{t,x.). 
From this it is clear that Eq. (|46)) can be expected to be 
valid only if S decreases with the increase in L, which 
forces the time of the development of the optimal fluctu- 
ation to coincide with L. In the opposite case (when S 
decreases with the decrease of L) there appears a possi- 
bility to decrease the action of the fluctuation we are con- 
sidering by making the time of its development smaller 
than L. Namely, if one makes in Eq. ([^^ a replacement 



L ^ j^L , A 7A 



(47) 



conserving relation ([jS)) . this leads to 5 ^ ^'^^'^S. There- 
fore, for d > 2 a consistent decrease in the size of the 
fluctuation and in the time of its development allow one 
to make S(F) arbitrarily small by choosing a sufficiently 
small 7. This suggests that the result can be ex- 
pected to be applicable only at < c? < 2, whereas at 
d > 2 the optimal fluctuation corresponding to the most 
distant part of the left tail has to be localized at small 
scales and its form has to be determined by the form of 
a cutoff. Without a cutoff the problem with d > 2 and 
(5-functional correlations is ill-defined. 

Note that at d > 2, the problem with (5-functional cor- 
relations is ill-defined also for another reason. Namely, at 
d > 2 the perturbative corrections to the viscosity v and 
other quantities acquire ultraviolet divergencies which at 
d < 2 are absent. Apparently, this is not a coincidence but 
another manifestation of the same phenomenon. There- 
fore, for d > 2 some ultraviolet cutoff must be introduced 
into the problem. One of the most natural ways to do it 
consist in assuming that the correlations of a random 
potential are characterized by a finite correlation radius. 



C. Finite-range correlations 

When random potential correlator U (x) [which we as- 
sume to be spherically symmetric, U (x) = J7(|x|)] is char- 
acterized by a finite correlation radius ^, the stationary 
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solution of Eqs. ([7]) and (fT4|) cannot be found exactly 
even at d = 1. However, it is clear from the form of Eq. 
(fT3| relating V and /i that when the soliton width A is 
much larger than ^, the actual solution has to be rather 
close to the solution for ^ = 0, the same being true also 
for the value of S{F). It follows from Eq. (f45|) that in 
terms of F the condition A 3> ^ corresponds to 



(48) 



It turns out that for the opposite relation between the 
parameters, ~F 3> F^, the stationary solution of Eqs. 
^ and (fT4|) also can be found rather accurately. As it 
is shown below, in such a case ji is localized in a region 
which is much narrower than ^, whereas both / and V 
change at the scales of the order of ^. In particular, it 
follows from Eq. that in such a situation the spacial 
dependence of the potential V^(x) just repeats that of 
(7(x), 

y(x) « -I7(x)e , (49) 
whereas the amplitude of V{'x) is determined by 



dx /i(x) 



(50) 



the overall strength of the negative potential source /^(x). 

For —F^F^ the viscous term in Eq. ([7]) can be ne- 
glected, which immediately gives that in the spherically 
symmetric stationary solution 



and 



df/dt w -U{0)e 



^] =2J[[/(0)-C/(r)]e, 



(51) 



(52) 



dr y/U{0) - U{r) . (53) 



where r — |x|, so that 
/(x) = /(O) 



In terms of the Schrodinger equation (|23|) the neglect of 
the viscous term in the stationary KPZ equation corre- 
sponds to nothing else but using the scmiclassical ap- 
proximation for the calculation of the ground-state wave 
function. 

Substitution of ^'(x) = cxp[— /(x)/T] with /(x) given 
by Eq. ([53]) into Eq. (f25|) demonstrates that at |x| <C ^, 



/i(x) c>c exp 



^2 1 



2A2 



(54) 



where A, the width of the region where the potential 
source //(x) is localized, is given by 



A{F) 



-,1/4 









1/4 



(55) 



When deriving this estimate we have replaced —Urr{0) 
by [/(0)/^2 and E by F/L. The result shows that the 
assumption A(F) ^ ^, which has been used above 
to obtain Eq. is indeed self-consistent as soon as 

Substitution of Eq. ((49|) into Eq. (fTS]) reduces the ex- 
pression for the action to a very simple form. 



S = 



2C/(0) 



(56) 



which is easily recognizable to those familiar with appli- 
cation of the optimal-fluctuation approach to a quantum- 
mechanical problem with finite-range correlations of a 
random potential"'^^ and after substitution of £' = F/L 
gives 



S{F) 



2U{0)L 



(57) 



The same temperature-independent answer can be also 
reproduced in terms of the Kardar-Zhang replica ap- 
proach (see Appendix [X)) . 

Thus we have demonstrated that for ^ > the most 
distant part of the left tail is Gaussian independently of 
the dimension. Since the width of the region where /i is 
localized grows with the decrease in —F, a crossover to 
some other regime must occur when this width becomes 
comparable with ^. In particular, for d < 2 and ^ <C 
the dependence of S* on F at —F ^ F^ has to be de- 
scribed by Eq. with a subsequent crossover to the 
universal regime discussed in Sec. IVIBI Naturally, the 
increase in ^ (or in d) leads to shrinking and subsequent 
vanishing of the region where S{F) can be described by 
Eq. ([46]) . On the other hand, when ^ is taken to zero F^ 
goes to infinity, which leads to the disappearance of the 
region with Gaussian behavior. 



D. A boundary from below 

It is worthwhile to emphasize that expression (|57p gives 
an exact boundary from below for the value of S{F) in 
the optimal fluctuation. This is so because the potential 
of the form 



(58) 



minimizes functional (|27p for the given value of y(x — 0), 
from where 



SiF)> 



1 



2U{0) 



dt[V{t,0)f 



(59) 



On the other hand, in a growing fluctuation of /(t,x) 
which has a spherically symmetric shape and an ex- 
tremum at x = 0, the absolute value of df{t,0)/dt is 
bounded from above by |l^(t,0)| because at the point 
of extremum the second term in the left-hand side of 
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the KPZ equation ([7|) vanishes, whereas the third term, 
—vV^f, has to have the same sign as df{t,0)/dt. This 
ahows one to conclude that 



S{F) 



> 



> 



1 



2U{0) 
2U{Q)L 



( dt 


\dfim] 


JO 





(60) 



Apparently, this inequality is reduced to equality only 
if (i) F(i,x) is of form ((58)) . (ii) the viscous term in the 
KPZ equation can be neglected, and (iii) df{t, 0)/dt does 
not depend on time. Since in the negative fluctuation of 
/ considered in Sec. lIVCl all these conditions are satisfied 
rather accurately, the action of this fluctuation is approx- 
imately equal to the boundary from below given by Eq. 

Note that the argument leading to the derivation of 
Eq. ((60|) is valid for both signs of F. Therefore inequality 
([60| has to be satisfied also in the far-right tail. 



V. FAR-RIGHT TAIL 

Our analysis has established that the optimal fluctu- 
ation corresponding to the left tail of Pl{F) has a very 
special shape which can be characterized by two different 
scales. Namely, the size of the area where the potential V 
is localized, A(F), is much smaller then the total size of 
the fluctuation A{F) - {-FL/J)^^"^, that is, the width 
of the area where / and u essentially deviate from zero. 
Apparently this property is closely related to the fact 
that inside a growing negative fluctuation of / the terms 
df/dt and (l/2J)(V/)2 in the functional. 



S{f} = 



1 

2U, 



Jo 



dl 
dt 



1 

2J 



{Vff - vV^f 



(61) 

defining the probability of a fluctuation in a system with 
a (5-correlated potential have to be of the opposite signs. 
This provides a possibility for their mutual compensa- 
tion in almost the whole volume of the fluctuation. It is 
clear that in the case of the right tail such a cancella- 
tion is impossible because in the substantial part of the 
optimal fluctuation df/dt has to be of the same sign as 
(1/2J)(V/)^. As a consequence, the optimal fluctuation 
corresponding to the right tail must have a shape which 
can be characterized by a single relevant length-scale, 
A+(F). 

This length scale can be estimated from the com- 
parison of df/dt ~ F/L with {l/2J){Vff - F^/JAl, 
which shows that A+ has to be of the same order as the 
total size of the optimal fluctuation with F < 0: 

1/2 



_(F)^A{-F)^[^ 



(62) 



Note that for (F) given by Eq. ([5^ the viscous term 
in the integrand of functional (pTjl can be neglected if F is 



large enough. This is precisely the reason why an estimate 
for A+ can be obtained by matching the two other terms 
in this integrand. A comparison of vV"^ f ^ vF/ A^ with 
(1/2J)(V/)^ ~ F'^/JA\ shows that the condition which 
allows one to neglect the viscous term can be written as 
F 3> 2Ji/ = T. Apparently this constraint is automati- 
cally fulfilled as soon as one considers the most distant 
part of the tail. 

Substitution of Eq. ([62| into the relation 



LA% 

~u7 



(63) 



following from the assumption that A+(i^) is the only 
relevant length-scale in the problem, gives then an esti- 
mate for the action determining the form of the far-right 
tail oiPL{F), 



S{F) 



p2+d/2 



(64) 



which naturally is independent of temperature. On a 
more formal level, the same relation can be obtained as 
a variational estimate from above. If one assumes, for 
example, that 



/(i,x) 



Ft 



■ exp 



2A2 



(65) 



and substitutes Eq. ([65|l into Eq. ([6111 , then for < d < 4 
the result of this substitution S'var(A+) (which for F ^ T 
is insensitive to the presence of the viscous term in the 
integrand) has a minimum with respect to the variational 
parameter A_|_. This minimum is situated at A^[F) sat- 
isfying relation ([S^ . whereas the value of 5'var[A+(F)] 
satisfies relation 

The important difference between the far-left and far- 
right tails is that in the far-right tail, the width of the 
region where the fluctuation of a random potential is lo- 
calized grows with the increase in \F\. In such a situation 
one can expect that the shape of the optimal fluctuation 
in the most distant part of the tail at finite ^ will be the 
same as for a (5-correlated potential. This requires the ful- 
fillment of the condition A+ » ^, that is, F ^ /F . 
Therefore, for a given ^ and sufficiently large L the region 
of the applicability of Eq. ([M|) will be extended to the 
whole non-universal part of the right tail. 

Although the minimum of 5'var(A+) with respect to 
A+ exists for any d in the interval < c? < 4, it follows 
from the form of Eq. (|64|) that this equation can be ex- 
pected to be directly applicable only at d < 2, exactly 
like in the case of the analogous expression for the far- 
left tail, Eq. (gni). For d > 2 Eq. dM]) (where L enters 
as the total time of the development of the fluctuation) 
predicts that the action can be decreased by making the 
time of the development of this fluctuation much smaller 
than L. According to Eq. (p^ this will be accompanied 
by the decrease in the size of the fluctuation. This sug- 
gests that at d > 2 the optimal fluctuation must have a 
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different structure, which has to be sensitive to the form 
of a random potential correlator at small lengths. 

If the first factor in the right-hand side of Eq. ([65]) is 
replaced by 

Fsinh(t/L+) 
sinh(L/_L+) ' 

which allows one to vary not only the characteristic size of 
a fluctuation A-|_ but also the time of its development L^, 
then for i7(x) cx exp(— and d > 2 the minimum 
of the action is achieved at A_|_ ~ ^ and L+ ^ J^"^ /F, 
which corresponds to 



S{F) 



UoJ 



(66) 



Note that at d = 2, the estimates given by Eqs. ((M)) and 
((66)) coincide with each other. Naturally, at the marginal 
dimension of d = 2 (where algebraic divergences are re- 
placed by logarithmic) some logarithmic factors may ap- 
pear in the expression for the action. 



VI. MODIFICATION OF TAILS BY THE 
RENORMALIZATION EFFECTS 

In terms of the Burgers equation parameters (the 
viscosity v ~ T/2J and the pumping force intensity 
D = Uo/2,P), Eq. can be rewritten as 



S{F) 



v^{^F/jf^ 

15 L^-dn 



(67) 



This estimate has been derived at d < 2 and ^ = and is 
applicable also at f > as soon as ^ ^ A. However, from 
the nature of the optimal-fluctuation approach it is clear 
that the range of the applicability of Eq. ([67]) is restricted 
also from the other side because in order to disregard the 
renormalization of any parameters by the nonlinearity 
the soliton has to be sufficiently narrow: A <C xq, where 
xq is defined by the relation 



„2-d 



15 



(68) 



At any d ^ 2 a;o is the only parameter with the dimension 
of X which can be constructed from T, J and Uq. In par- 
ticular, in the case of d < 2 and small ^ we are discussing 
now, xq is the length scale at which the perturbative cor- 
rections to v and D become comparable with the bare 
values of these parameters. 

Thus, at A ^ xo the renormalization effects become 
important. In such a regime the probability of a large 
negative fluctuation of F is determined not by a single 
fluctuation (and small deviations from it) but by a rel- 
atively wide class of fluctuations, the summation over 
which can be taken into account by analyzing an optimal 
fluctuation in a system with renormalized parameters. 



Since in all the cases we consider the optimal fluctua- 
tions are quasi-stationary (see below) and well localized 
at a particular length scale, this can be done by replac- 
ing all parameters in Eq. (|67p by their effective values at 
the corresponding length scale and zero frequency.^ How- 
ever, it is well known that only v and D are subject to 
renormalization, whereas the amplitude of the nonlinear 
term in the KPZ equation ([7]) (and, therefore, the coeffi- 
cient J) cannot be renormalized as a consequence of the 
Galilean invarianccpi^ 

From the continuity it is clear that when the instanton 
is not too narrow, the approach relying on using Eq. ([^7]) 
with renormalized parameters can be also expected to 
work even at d > 2 [where Eq. ([S7|) has no region of 
the direct applicability] as soon as the parameters of the 
system correspond to the same phase as at d < 2 [namely, 
the strong-coupling phase in which the fluctuations of 
f{t, x) in a stationary situation are divergent, see Eq. 
below]. At d > 2 this requires to have xo/£, > K{d), that 
is, the temperature T should be lower then some critical 
value Tc{d)r^ which tends to infinity when d 2 -I- 0. 
In the weak-coupling phase, that is at T > Tc{d), typical 
fluctuations of f{t,x.) in the stationary situation can be 
described by neglecting the non-linear term in the KPZ 
equation ([7]) . However, the form of the most distant parts 
of the tails of (F) is insensitive to the relation between 
T and Tc(d) and in both phases has to be given by Eqs. 
dSZl) and (ESI). 

To describe how the renormalization effects change the 
shape of the tails of Pl (F) in the regime when they are 
important (which corresponds to the universal parts of 
the tails in the strong-coupling phase), we flrst have to 
review some known properties of the stationary solution 
of the KPZ model in the strong-coupling regime. 



A. Stationary solution of the KPZ model 

In a stationary situation the divergence of fluctuations 
in the strong-coupling phase of a KPZ system is alge- 
braic. Their behavior at large scales in space-time can be 
described by two fundamental exponents, ^^'^^ 



([/(i,x)-/(t',xT)(x|x-xf'^5 



\t-t'\ 



(69) 



Here x = xW is the roughening exponent character- 
izing the equal-time interface fluctuations, z = z{d) is 
the dynamic exponent describing the scaling of the re- 
laxation time with the length-scale, whereas the function 
g{a) has a finite limit at a and diverges as a'^^^^ 
when a — > oo. It is well known^* that the existence of 
the Galilean invariance imposes 



z + x = 2. 



(70) 



At d = 1 the value of the exponent x = 1/2 is known 
exactly because the equal-time correlator of /(t,x) in a 
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system with (5-functional correlations of a random po- 
tential has to be exactly the same as in the absence of 
the non-linearity.— This property is a consequence of the 
fluctuation-dissipation theorem, which is obeyed by Eq. 
([7]) only at d = 1. At d 7^ 1 the values of the exponents 
z and X ^'^^ known only from approximate or numerical 
calculations. In terms of the directed polymer problem 
the dependence corresponds to 



([x(^)-x(^')n(x(^-^')'/^ 



(71) 



which shows that = Xjz plays the role of the rough- 
ening exponent for the transverse displacements inside 
an infinite polymer and therefore cannot be smaller than 
l/2,i2^ from where z <1. 

A natural way to describe the effective renormalization 
of V and D by the nonlinearity consist in introducing^'^ a 
generalized viscosity i/(ci;, q) and a generalized pumping 
intensity -D(a;,q) defined by the relations 

G(c^,q) = {-iu^v(u:,^)qW (72) 
C(c^,q) = 2j2|G(c^,q)pI?(c^,q), (73) 

where G(io^c\) and C(w,q) are, respectively, the Fourier 
transforms of the response function and of the two-point 
correlation function of /(t, x). The form of Eqs. (|72|) and 
(175)) corresponds to the replacement of the considered 
non-linear system by a linear system with the same form 
of G(ti;,q) and G(ti;,q). 

The compatibility with the behavior described by Eq. 
requires then that at small enough g. 



lim v{lo^ q) oc q 



lim -D(w, q) oc q 



This suggests that the behavior of low-frequency fluctu- 
ations with typical or smaller amplitude can be quali- 
tatively described by using an effective viscosity Vce{R) 
and an effective pumping intensity DesiR) which alge- 
braically depend on a length scale R, 



VesiR) 



R 



2-z 



D,s{R) ^ D 



R 

an 



4+d-3z 



(74) 



where in accordance with Eq. ([70| we have replaced x by 
2 — 2. As a convenient way of describing the amplitudes 
of i'cff(-R) and Dce{R), we have introduced in Eq. ([7i|) 
two new length scales, a,^ and an- For d = 1 and ^ ^ xq 
both and ao can be expected to be of the order of 
Xq, because in such a situation xq is the only relevant 
length in the problem. However for ^ 3> and/or d > 1 
these two length-scales do not have to be of the same 
order. Since both ly and D increase under the renormal- 
ization, Eqs. ([74]) can be expected to be applicable only 
for R^ a^^ajj. 

In scaling regime, when Ves{R) and Deff(-R) behave 
themselves in accordance with Eqs. ([74| . both these 
quantities have no direct relation to their bare values, 
1/ and D. Their origin can be traced to the effect of fluc- 
tuations with shorter wave lengths than the given length 



scale R. In particular, it follows from the structure or 
the KPZ equation ([7]) that at the length-scale R the role 
of the effective random potential is played by the devia- 
tion of — (J/2) (u^);? from its average value, — {J / 2) (u^) , 
where (. . denotes spatial averaging over a region with 
a linear size of the order of R. From this the value of 
Dcs{R) can be estimated as 



DesiR) 



dr 



\r\<R 



dr 



u''{t,x.)u''{t + T,x + r) 



fcsiR) 



(75) 



In Eq. (|75p we have assumed that the integration over 
dr can be replaced by the multiplication by the factor 
~ T(r), where r(r) r^/uesir) the characteristic re- 
laxation time which can be associated with the length 
scale r, whereas the result of the integration over dr has 
been estimated assuming that as a consequence of the 
universality for any length scale there exist only one char- 
acteristic velocity scale which can be associated with this 
length scale (in other terms, there is no anomalous scal- 
ing). We have chosen as such a velocity scale the typical 
velocity, utyp(i?), defined by the relation 



(i?) 



"typ 

Substitution of Eq. 
relation 



u 



DcsjR) 
v.s{R)R'' 



(76) 



76)1 into Eq. (|75)) then gives the 



DesiR) 



R 



2-d 



(77) 



whose structure is analogous to that of Eq. ([68]). The 
consistency between Eqs. (|77)) and (|70| confirms the cor- 
rectness of assumptions which have been used for the 
derivation of Eq. (|77)) . In terms of the length scales Qi, 
and aD introduced above, see Eqs. (|7^ . relation ()77p can 
be rewritten as 



Xq 



an , 



which for d = I (when z = 3/2) is reduced to 



(78) 



(79) 



When the dynamics of fiuctuations at i? ~ ^ is dom- 
inated by wave breaking, the value of Mtyp(C) can be 
estimated as a characteristic velocity ^ {Dt^/S^^Y^"^ , 
which is created by a random force with character- 
istic length-scale ^ during the time ~ C/^^ re- 
quired for breaking of such a fluctuation, which gives 
^ (D/S^^^'^y^'^ . A comparison of this estimate with 
Eqs. (|75)) and suggests that in such a regime 

DcsiO ^ -Oj that is, od ~ ^. At d < 2 we expect this con- 
clusion to be applicable when ^ xq, whereas at d > 2 - 
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in the whole region of the existence of the strong-coupHng 
phase. 

It follows from the definition of utyp(-R) that with the 
increase in R the value of utyp(i?) has to decrease. A 
comparison of Eq. ((76)) with Eqs. ([74)) allows one then to 
conclude that z has to be larger than 1. 



B. Universal part of the left tail 

After replacing in Eq. (|45| T/2J = ly by i'off(A), one 
obtains a relation which allows one to find that in the 
regime when the renormalization effects are important 
the estimate for the soliton width A acquires a form 



A = A(F) 



-Fal) 



(80) 



A substantial change of A in comparison with what is 
given by Eq. (|45p means that in the regime we consider 
now the probability of a large negative fluctuation of F 
is determined not by a narrow vicinity of the fluctuation 
which minimizes the original action (like it happens in 
the more distant part of a tail), but by a wide vicinity 
of an essentially different fluctuation whose dominance 
is ensured by a factor related to the integration over its 
vicinity. In the framework of a renormalization group ap- 
proach this factor is effectively taken into account when 
one is replacing different parameters by their renormal- 
ized values. 

An estimate for the action can be then obtained by 
making in Eq. ()67[) a replacement 



(81) 



1^ ^ i^eff(A) , D ^ DMA) 



with A given by relation With the help of Eq. (|77)) 
the result of this substitution can be reduced to the form 



with exponent 



SiF) - { — 



11 = 1- 



2(2-1) 



(82) 



(83) 



which depends on d only through the dynamic exponent 
z = z{d) but not explicitly. Here 



(84) 



plays the role of a characteristic free-energy scale whose 
dependence on L is described by the exponent 



1 2 

w = 1 = - - 1 

z 



(85) 



The universality hypothesis for the directed polymer 
problem?^ (or, more generally, for the collective pinning 



problem-^) suggests that has to be of the same or- 
der as a characteristic elastic energy Ec\ ^ J{Sx)'^/L, 
where the dependence of the characteristic transver- 
sal displacement between the two ends of a polymer, 
Sx = |x(i = L) — x(i = 0)1 (X L'', on its total length L is 
described by the roughening exponent C so that 



= 2C - 1 



(86) 



A comparison of Eq. ([86)1 with Eq. ([85)) demonstrates 
that the fluctuations of dx are described by the same 
roughening exponent C = 1/z as fluctuations inside an 
infinite polymer, see Eq. I|7ip . in full agreement with what 
one expects from the universality. This consistency can be 
considered as an additional confirmation of the validity of 
the set of assumptions which have been used for obtaining 
Eq. 

Note that the list of these assumptions includes the 
conjecture that the system evolves sufficiently slow, so 
that at relevant length scales it can be considered as al- 
ready equilibrated, which is a necessary condition for us- 
ing Eqs. ([74)) . For this the total evolution time L has to 
be much larger then the characteristic relaxation time 
t(A) ~ A^/i'eff(A) which can be associated with the 
length scale Since in terms of A{F) and L relation 
(l82l) can be rewritten as 



SiF) 



t^cff(A) 
A2 



L 



L 



r(A) 



(87) 



the constraint L ^ t{A) is equivalent to S{F) ^ 1 and, 
accordingly, is automatically fulfilled as soon as one is 
dealing with the tail. 

It is also important that the effective viscosity and ef- 
fective pumping intensity given by Eqs. (|74p and follow- 
ing from the form of the correlation function (|69p can be 
used for the description only of typical (or more weak) 
fluctuations. The comparison of the characteristic veloc- 
ity of the flow created around the instanton. 



Up 



JL 



1/2 



with utyp(A), the typical velocity of equilibrium fluctua- 
tions at the length scale A [see Eq. ([76)) ]. demonstrates 
that in the considered case both quantities are of the 
same order, and therefore, the approach based on using 
Eqs. ([74)) with i? ^ A is indeed justified. This allows us to 
conclude that our instanton is created by fiuctuations of 
the effective random potential whose amplitude is typical 
for their length scale. In such a situation the only reason 
why the probability of the instanton is small is that the 
signs of these typical fluctuations have to be same in all 
L/t{A) independent time intervals of the length t(A). 
This provides a qualitative explanation why the expres- 
sion for the action can be reduced to a very simple form 
5-i/r(A). 

At large values of —F the range of the applicability of 
Eq. ([5^ is restricted by the constraint A{F) ^ a^,au^ 
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which is required for making replacement (jSip . In partic- 
ular, when d < 2 and ^ ^ (so that ^ ajj ~ xg) one 
can expect that at A(F) ~ xq, that is, at 



Jxq 



(89) 



a crossover takes place from dependence ([5^ to depen- 
dence pS)) . 

On the other hand, in situations when ^ (or d) is too 
large for dependence to have any range of applicabil- 
ity, one could expect to have a direct crossover between 
dependences (|57p and (|82p. However, the range of the 
applicability of Eq. ((57)) describing the far-left tail corre- 
sponds to A(F) <^ min[a{,a£)] and of Eq. ([5^ describing 
the universal regime to A{F) ^ max[a^,a£i]. Since we 
expect that in a general situation the two length scales, 
a J and ao, are essentially different, we have to admit 
the existence in such a case of an intermediate region in 
F, where the form of the left tail of Pl{F) cannot be 
established without further investigation. 

As it has been already mentioned above, at d = 1 the 
value of the exponent z is known exactly. Substitution of 
z = 3/2 and ^ xq into Eq. (EH) and Eq. ^ then 
reproduces an estimate for S{F) which up to unknown 
numerical factor coincides with Eq. (|55|) for S{F) in the 
far-left tail. This shows that for d = 1 and ^ ^ a;o the 
dependence of S on all parameters in the universal part of 
the left tail is exactly the same as in its non- universal part 
at Fc <C —F <C F^. In this particular case at —F ^ Fc 
only a numerical coefficient in the dependence (|82p can 
experience a crossover. 

For d — 1 and f xq, substitution of Eq. ((79|) with 
an ^ £, into Eq. ([M)) reproduces an estimate for F^,{L) 
which has been obtained by Nattermann and Renz'^^ from 
scaling arguments complemented by the assumption that 
at low enough temperatures F*(L) has to be temperature 
independent, and follows also from the replica-symmetry- 
breaking analysis of Ref. 28. 

For d ^ \ the value of the exponent rj in the universal 
part of the left tail, rj- = z/[2{z — 1)], does not coincide 
with its value in the far-left tail, where it is given by 
2 — d/2 [see Eq. (|^5))]. Note that expression ([5^ decreases 
with the increase in L as long as 7y_ > 1, that is, 1 < 
z < 2. This means that in the universal part of the left 
tail the condition which is necessary for the possibility of 
having a macroscopic optimal fluctuation (whose size is 
much larger then is changed from (i<2tol<z<2. 
On the other hand, when the renormalization effects are 
taken into account, the condition < d < 4 required 
for having a minimum of S with respect to A (see Sec. 
HVB]) is replaced by 4/3 < z < 2. Thus, the range of the 
applicability of Eq. ([5^ is not restricted to < c? < 2 
(as in the case of the analogous expression for the far-left 
tail) but extends itself to the whole region of parameters 
where the strong-coupling phase does exist and z > 4/3 
(the condition z < 2 always has to be fulfilled, see Sec. 
I VI Ap . Note that for d = 1 the value of C = 1/z is equal 
to 2/3 and according to numerical simulations goes down 



with the increase in d.^ This means that the restriction 
z > 4/3 is fulfilled for any physical dimension. 

C. Universal part of the right tail 

One could expect the approach based on the applica- 
tion of the replacements (|5T|) to be applicable also for the 
description of the universal part of the right tail. How- 
ever, it turns out that in this case the situation is more 
complex. This can be understood by comparing the size 
of the optimal fluctuation A+(F), given by Eq. ((62)) . with 
the length scale i?*(i^) at which the typical velocity of 
equilibrium fluctuations wtyp(i?*), given by Eq. ([75)) . be- 
comes comparable with 



Up 



F 

7Ki 



F 

7l 



1/2 



(90) 



the characteristic velocity inside the optimal fluctuation 
with the size A+(F). In Eq. (j90p we have used the es- 
timate for A+{F) given by Eq. (p^ . which has led to 
exactly the same estimate for wj? in terms of \F\ as in 
the left tail [see Eq. (j88p ] . This means that in both tails 
R^,{F) has to be of the same order. On the other hand, 
in Sec. IVIBI we have established that in the left tail the 
relation up ^ Utyp(-R*) holds precisely when i?* ~ A(i^). 
This allows one to conclude that in the right tail. 



K{F)^A{-F), 



(91) 



where A(F) is the instanton width in the left tail given 
by Eq. dHOl)- 

Accordingly, for the creation of the optimal fluctua- 
tion whose size A+{F) is much larger than R^,{F) (as 
it is required in the case of the right tail), the fluctua- 
tions of the effective random potential with length-scale 
A^{F) should have amplitudes much larger then typical. 
Naturally, the probability of such fluctuations is strongly 
suppressed and cannot be estimated by using Eqs. ((8T|) . 

The most effective way of formation of a fluctuation 
whose amplitude up substantially exceeds the typical ve- 
locity of fluctuations at the corresponding length-scale 
consists in formation of a set of fluctuations with smaller 
length scales, such that for them the amplitudes of the 
order of up are typical. This means that the length-scales 
of these fluctuations should be of the order oiR^,{F), and, 
accordingly, the estimate for the action should include an 
additional factor (A^/i?*)"* which takes into account the 
need for the spatial coherence of these fluctuations. This 
leads to 



S{F) 



L 



r(i?*) 



v+ 



(92) 



where F^, is the same characteristic free-energy scale as in 
the universal part of the left tail [see Eq. ([84]) ] , whereas 
exponent 77+ is given by 



{l + d)z 
2(^-1) 



(93) 
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In terms of the renormalization approach exactly the 
same result is obtained if the renormalization is stopped 
not at the scale A_|.(F), corresponding to the total size 
of the optimal fluctuation, but at a smaller scale i?* (at 
which the fluctuations stop to be strong enough for in- 
ducing the renormalization), that is, by using Eq. ([64]) 
with the replacement 



4+d-3z 



(94) 



where ~ A{~F). Since the value of Deft(^*) does 
not depend on A+, the condition for the existence of 
a minimum of S with respect to A_|_ remains the same 
as has been found when deriving Eq. ((64)) . < d < 4. 
On the other hand, in the universal part of the right 
tail the condition required for the possibility of having 
a macroscopic optimal fluctuation (whose size is much 
larger then ^) is changed from d<2tol<z<2, which 
in the strong-coupling phase anyway has to be fulfilled 
[see Sec. IVI A] . Therefore, the range of the applicability 
of Eq. ([g^ is restricted from above not by d = 2 (as in 
the case of the analogous expression for the far-right tail) 
but by d = 4. 

Note that in contrast to exponent rj- given by Eq. 
exponent 77+ depends both on z and d. However, the ratio 
of these two exponents does not depend on z, 



^ = l + d. 
V- 



(95) 



and therefore is known exactly. The fact that in the 
regime where the renormalization effects are important 
both tails of the free energy distribution function incorpo- 
rate the same characteristic free-energy scale Fs, confirms 
that this regime corresponds to studying the universal 
form of this distribution function. 

A comparison of Eq. ((62)) with Eq. ([80| allows one 
to verify that the condition A_|.(F) ^ R^,{F), on which 
we have relied when deriving Eq. (j92p . is equivalent to 
S{F) ^ 1, and therefore is always satisfied as soon as we 
are dealing with the tail. Another condition whose fulfill- 
ment is required to justify replacement (|94p is related to 
the quasistationarity of the problem. Namely, the total 
evolution time L has to be much larger than the charac- 
teristic relaxation time r(i?*) ~ Rl/vcs{R*) which can 
associated with the length scale R^{F). For i?*(F) ^ 
A{—F), this condition is also reduced to S{F) 1. 

From the side of large F the range of the applicabil- 
ity of Eq. ([^ is restricted by the constraint i?* a^i , 
whose fulfillment is also required for making replace- 
ment In particular, at d < 2 and £, ^ xq (when 
o-D ^ xq), the crossover between dependences (j92l) and 
([M)) can be expected to occur at F ^ Fc, where Fc is 
given by the same expression [Eq. ([55)1 ] as in the left tail. 
On the other hand, at d > 2 the crossover between de- 
pendences (|92)) and ((66|) has to take place while i?* (F) is 
still much larger than ^. In this situation we expect that 
the two contributions to Pl{F) [one from the "macro- 
scopic" instanton, corresponding to dependence (|92p and 



the other from the "microscopic" instanton correspond- 
ing to dependence can coexist with each other and 
the crossover has to occur when they become comparable 
with each other. 



VII. CONCLUSION 

In the present work we have studied the form of the 
tails of the free-energy distribution function Pl {F) in the 
directed polymer problem both for a (5-corrclatcd random 
potential and for the case of a finite correlation length ^. 
In all regimes that we have investigated the tails have a 
stretched-exponential form. 



lnPi(±F) 



F*(L) 



(96) 



with F^,{L) (X L^^ and therefore can be characterized by 
the two exponents whose values depend on the dimen- 
sionality of the space in which the polymer is imbedded. 
We use letter d to denote the transverse dimensionality 
of this space, that is, the number of components of the 
displacement vector u. 

For sufficiently large fluctuations of F the form of the 
tails of Pl{F) is determined by the form of the most op- 
timal fluctuation of a random potential which is sufficient 
for achieving a given value of F. For a (5-correlated ran- 
dom potential and d < 2 the minimization of the action 
corresponding to such a fluctuation allows one to show 
that in the far-left tail 



4-d 



2- d 



V- = 



(97) 



The same values of tj- and w_ have been obtained by 
Monthus and Garel^^ by constructing a generalization 
of the Imry-Ma scaling argument (based on a disorder- 
dependent Gaussian variational approach introduced in 
Ref. l30h . However, the approach of Ref. [2^ leaves one in 
doubt on what is the range of its applicability (and if 
such a range exists at all), whereas the methods used in 
this work allowed us to establish that the exponents ([^7]) 
are applicable in the most distant part of the left tail 
corresponding to the nonuniversal regime. 

At d > 2 the problem with strictly (5-functional corre- 
lations of a random potential becomes ill-deflned, so it 
becomes necessary to introduce some regularization. The 
natural way of doing it consists in assuming that a ran- 
dom potential correlations are characterized by a finite 
correlation radius In the case of ^ > one finds that 
in the most distant part of the left tail the size of the 
optimal fluctuation of a random potential has to be com- 
parable with ^ and the values of the exponents become 
superuniversal, that is, not dependent on d. 



^- =2, 



= 1/2, 



(98) 



For d < 2 and not too large ^ one can expect to have a 
crossover from regime (|M)) to regime (jHT])- 
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The application of the optimal-fluctuation approach to 
the analysis of the right tail shows that for d < 2 the most 
distant part of this tail is described by 



2-d 



(99) 



In contrast to the case of the left tail, the form of the most 
distant part of the right tail is insensitive to whether ^ 
is zero or finite. On the other hand, for d > 2 the size of 
the optimal fluctuation again becomes determined by ^, 
which leads to the change of the exponents to 



LU:=0. 



(100) 



Note that the value of uj+ given by Eq. (f99|) corre- 
sponds to the value of the roughening exponent. 



Cf 



(101) 



which is known as "Flory exponent"— and follows from 
simple scaling arguments of Refs. [U as well as from the 
Gaussian variational calculation of Mezard and Parisi^- 
incorporating a hierarchical replica-symmetry breaking. 
Our analysis has revealed that this scaling analysis 
(which insofar has been assumed to be of little relevance, 
because it cannot reproduce the exactly known value of 
C = 2/3 at d = 1) in reality is applicable for the descrip- 
tion of the most distant (non-universal) part of the right 
tail of Pl{F). However, it still remains unclear whether 
the appearance of the same exponent in the variational 
calculation of Ref. [S^ (based on the maximization of the 
variational free energy of a system with i = oo) is not 
more than a coincidence. 

If the parameters of the system correspond to the 
strong coupling phase, the decrease in \F\ makes the 
optimal-fluctuation approach no longer directly applica- 
ble because the size of the optimal fluctuation becomes 
too large (or its amplitude becomes too small) to neglect 
the renormalization of the parameters of the system by 
fluctuations. In such a situation, a consistent inclusion 
of the renormalization effects into account allows one to 
express the exponents in terms of the roughening expo- 
nent C = ^/z describing the behavior of displacement 
fluctuations inside an infinite polymer [see Eq. (|7ip ]. For 
universal parts of left and right tails, one obtains, respec- 
tively, 
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1 



2(1-0 



2(1-0 



(102) 



Not unexpectedly, one finds that the value of lo is the 
same for both tails and is equal to 2C — 1, as it could 
be expected from the universality. Quite remarkably, the 
ratio 77+/?7_ = 1 + d does not depend on 

The value of C, is known exactly only at d = 1, where 
C = 2/3. In this case the values of ?]_ — 3/2 and 
77+ = 3 which follow from Eqs. (|102p are in perfect 
agreement with the exact solutionil of the polynuclear 



growth (PNG) model, which is accepted to belong to the 
same universality class. In terms of the directed poly- 
mer problem the PNG model corresponds to the Pois- 
son distribution of identical pointlike impurities and a 
rather peculiar limit of vanishing elasticity, J — Q, and 
zero temperature] ^^'^'^ For this model the form of the 
distribution function Pl{F) in the universal regime, as 
well as the scaling function g{a) entering Eq. (j69p . is 
known exactly j^^i'^^ The consistency between our results 
and that of Ref. [ll| confirms that the directed polymer 
problem defined by Eq. ^ and the PNG model indeed 
belong to the same universality class. 

The nonuniversal tails in the PNG model have been an- 
alyzed in Ref. [ssl . Naturally, in the nonuniversal regime 
even the models belonging to the same universality class 
can have different tails. The difference is especially evi- 
dent in the case of what we call the far-right tail because 
in the PNG model the energy is by definition bounded 
from above and therefore its distribution function has 
to vanish for large enough positive fluctuations. On the 
other hand, it follows from Ref. [ssl that in the PNG model 
the far-left tail is described by S{F) oc F\n{—F/L) and, 
thus, also has nothing in common with the far-left tail of 
the model considered in this work. 

In terms of the exponent a; = 2^ — 1 Eqs. (|102p can be 
rewritten as 



V- 



1 



1 



(103) 



Our results demonstrate that in model ^ the analogous 
relations are fulflUed also in non-universal regimes (where 
uj is not obliged to coincide with 2^ — 1 and be the same 
in both tails) as soon as the size of the optimal fluctua- 
tion is comparable with the total length of a string. For 
the far-left tail this has been known-'-^^ from the Kardar- 
Zhang replica approach. Recently both relations (|103p 
have been derived by Mothus and Garel"^^ with the help 
of a recursive procedure for the zero-temperature prob- 
lem on a hierarchical diamond lattice whose effective di- 
mension is equal to d. These authors have also suggested 
that the same relations can be expected to hold on all hy- 
percubic lattices. Although, in our opinion, the argument 
accompanying this proposal does not take into account 
some important differences between hypercubic and hier- 
archical lattices, the results derived in this work conflrm 
its correctness. 
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APPENDIX A: THE REPLICA APPROACH 

The replica approach to the directed polymer problem 
is based on calculating the moments Z„ = of the 
distribution of the partition function Z = z(L, 0) and 
allows one to find the far-left tail of Pl{F) without re- 
lying on the analytical continuation of n to 0. Kardar'* 
was the first to notice that for any integer n > 1 (and 
large enough polymer length L) Zn with an exponential 
accuracy can be approximated as 

Z„ « eM-Eo{n)L/T] , (Al) 

where Eq (n) is the ground state energy of the quantum- 
mechanical Hamiltonian, 

„2 " n n 

^" = - 2J E - 2T ^ ^ ^^^-^ 

a— 1 a— 1 6^1 

describing n bosons whose mass is equal to J (with T 
playing the role of h) and interaction to — J7(x)/r. In 
a (1-1-1) -dimensional system with a ^-correlated ran- 
dom potential, U{x) = UqS{x), the ground-state wave- 
function for the Hamiltonian (IA2p and its energy, 

i?o(n) = -^n-^n(n^-l), (A3) 

can be found exactly*^ This gives 

Z^cxexpl^^n^L^ , (A4) 

where the linear in n term in EQ{n) has been omitted, 
because it can be eliminated by a constant shift of the 
potential V{t,x) in Eq. 0. Note that the form of Eq. 
(|A4[) does not depend on the initial condition. The choice 
of the initial condition manifests itself only in the form 
of a prefactor which in the first approximation can be 
ignored. 

Since Z" = exp(— nF/T), Z„ can also be expressed 
in terms of Pl{F), the distribution function of the free 
energy F= /(L,0) = -TlnZ: 

/ + 00 
dFPi(F) exp(-n^^/T). (A5) 

It is clear that the integral in Eq. (|A5p can be convergent 
for any integer n> 1 only if the left tail of Pl (F) decays 
faster than exponentially. In such a case the integration 
in Eq. (IA5P can be performed with the help of the saddle 
point approximation. A comparison of the two expres- 
sions for Zn [Eq. (|X4)) and Eq. (|X5)) ] has led Zhang to 
conclusion^ that the algebraic growth of In Z„ cx Ln^/" 
at large n can be recovered only if one assumes that 
when —F is large S{F) = — ln[Pi(F)] is proportional 
to {-F/L'^)'^ with 77 = 1/(1 - uj).^ Thus, 1/w = 3 corre- 
sponds to 77 = 3/2. 

It is not hard to check that expression ([38|) for S{F) 
derived in Sec. lIVAl allows one to reproduce not only the 



power of n but also the full expression (jA4p . including 
the coefficient in front of rr" . For such a form of the tail 
the integral in Eq. (jA5| at positive n is dominated by the 
vicinity of the saddle-point situated at 

1TJ2t 

F = -F,{L)n\ F,{L)^^, (A6) 

where the full expression standing in the exponent, 

S{F) = ^S{F) - nF/T , (A7) 

has a maximum with respect to F. Substitution of Eq. 
(K6\\ into Eq. (K7\ leads then to Eq. (K^ . whereas cal- 
culation of d'^S{F)/dF'^ = l/{2TFcn) allows one to ver- 
ify that the condition for the applicability of the saddle- 
point approximation has a form {Fc/T)n^ ^ 1 and is 
automatically fulfilled for any integer n > 1 as soon as 
Fc ^ T , which anyway is required for the applicability 
of Eq. jSH). 

Since the explicit expression for Z„ given by Eq. (|A4[) 
can be derived only at integer n > 1, the region of 
the applicability of the replica approach is restricted to 
— F 3> Fc{L), that is, coincides with the region where 
the left tail can be found with the help of the optimal- 
fluctuation approach (see Sec. IIV Ap without taking into 
account the renormalization effects. However, the anal- 
ysis of Sec. IVIBI has demonstrated that in the (1 + 1)- 
dimensional systems with (5-functional correlations the 
form of the left tail in the region F^, — F Fc, where 
the renormalization effects have to be taken into account, 
remains qualitatively the same as for —F ^ Fc- In terms 
of the replica approach this means that the analytical 
continuation of the partition function of model (jA2p even 
at n < 1 behaves itself as if it was dominated by the con- 
tribution from the same state as at n > 1, although at 
71 < 1 this state no longer has the lowest energyi^ The 
reasons for that still remain to be elucidated. 

For ^ > the form of the ground state of the Hamil- 
tonian (jA2p cannot be found exactly even when d = 1. 
However, for 1^ ^ a;o the value of Eo{n) in a (1 + 1)- 
dimensional system with ^ > can be found rather ac- 
curately for any integer n > I, because in this regime the 
typical distance between bosons is much smaller then the 
radius of their interaction^^ and thus the main contribu- 
tion to Eo{n) is given simply by —U{0)n^/2T. A com- 
parison of Z„ - exp[f7(0)L7iV2T2] with Eq. ([All then 
immediately leads to the conclusion that the form of the 
far-left tail must be described by = F^ /2U{0)L, as it 
has been already derived in Sec. IIVCI with the help of 
the optimal-fluctuation approach. From the origin of this 
result it is clear that it has to be applicable for the de- 
scription of the most distant part of the left tail at finite 
^ for any d. 

The two contradicting attempts of generalizing the 
replica approach to d > 1 in the regime when the main 
contribution to Eq^u) is determined by the full form of 
[/(x) (and therefore cannot be found without introducing 
some additional approximations) have been undertaken 
by Zhang'^^ and Kolomeiskyi^ 
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APPENDIX B: FIXED INITIAL CONDITION 

If at i = the polymer is fastened at a; = 0, the ini- 
tial condition for the partition function z{t,x) has to be 
written as 



z(0, x) = XdS{x) , 



(Bl) 



where a dimensional factor Xd (with the same dimension- 
ality as x) has to be inserted to make z{t, x) a dimension- 
less quantity. This factor has no physical meaning and 
must drop out from all physical quantities. 

In the absence of pumping the initial condition (|Bip 
leads to 



z{t, x) = 



(47ri/i)i/2 



exp 



X 



= 7S°\t,x) (B2) 



and 



Jx^ 
~2F 



W-^^f^"\t,x), (B3) 

Z Xi 



which in terms of tt = V xfit, x)/ J corresponds to 

u{t,x) ^ x/t = u'^"\t,x) . (B4) 

Accordingly, the initial condition for f(t, x) can be for- 
mulated as 



lim 



fit,x)-f<^°Ht,x) 



0. 



(B5) 



In such a situation it seems to be convenient to count off 
the free energy from its value in the absence of pumping 
introducing 



fit,x)^fit,x)-f^'Ht,x), 



(B6) 



(i) because f{x,t) in contrast to f{t,x) does not depend 
on Xd and (ii) because this allows to rewrite the initial 
condition (IB5I) as 



f{0,x)=0. 



(B7) 



Accordingly, in the case of fixed initial condition F should 
be redefined as 



^^ = /(L,0) = /(L,0)-/(")(L,0) 



(B8) 



It is clear that the stationary solution ([5^ in no way 
resembles the initial condition (|B5[) . However, like in the 
case of the free initial condition, it turns out possible to 
modify this solution without increasing the action in a 
way which eliminates the inconsistency with the initial 
condition. 

This solution has to interpolate between the soliton 
at small x and the dependence f{t, x) — (i, x), which 



has to survive in the regions that in a given moment 
are too far from the soliton to feel its presence. One can 
expect that the time of the formation of a solution, ti, 
will be of the same order as for free initial condition, ti ~ 
i^/v'^, because a.tt ^ ti and a; ~ A the velocity of the flow 
(IB4p produced by fixed initial condition is much smaller 
then the velocity v ^ v/ A produced by the soliton with 
the width A. Then a.t t ^ ti and A <^ x <^ vt f has to 
be of the form 



.fit,x) = ~ — it-t,) + Jv\x\ + -ln—^ 



f'^'\t,x) 



(B9) 

where the constant has been chosen in such a way that 
the extrapolation to t — ti gives /(ti,0) = /^°-'(ti,0). 

A comparison of Eq. (|B3p with Eq. (|B9p shows that 
the crossover between these two dependences has to take 
place in the vicinity of |a;| — vt. Since we assume that at 
|x| A the potential is absent, this crossover has to be 
described by 



fit,x) = f^'Ht,x) ~T\nziit,\x\ - vt) , 



(BIO) 



where the function zi{t,x) is the exact solution of the 
diffusion equation zt = vz^x which at i = t\ = 
^1 exp[(Jw^/2T)ti] ~ ti [that is, when /(''^(t, wi) coincides 
with /(°)(t,tii)] smoothly interpolates between zi « 1 at 
-x > A and zi « exp[- Ja::^/2Tt'i] at a; > A. There- 
fore, the asymptotic behavior of zi{t^ x) at large t can be 
found by assuming 

zi{t'^,x) ^ 0{-x) + e{x) exp[~Jx^/2Tt[] , (Bll) 

where 9{x) = ■^[l+sign(a;)] is the step-like function. How- 
ever, for establishing the relation between S and F the 
exact form of zi{t,x) is irrelevant. We only have to be 
sure that at x ~ —vt this quantity is very close to 1, and 
this is satisfied as soon a.s t ^ t[ . 

The main difference with the case of free initial condi- 
tion appears in the relation between A and F. Subtrac- 
tion of Eq. (|B3p from Eq. (jB9p shows that for fixed initial 
condition Eq. (|37p should be replaced by 



T2 T L 

F = —L — 

2JA2 2 t. 



(B12) 



However, for L ^ ti the second term in the right-hand 
side of Eq. ([57]) is much smaller then the first one (which 
is of the order oiTL/ti) and therefore can be neglected. 

This allows one to conclude that in the case of fixed 
initial condition the main contribution to the action has 
exactly the same form as for free initial condition. The 
derivation above can be easily generalized for the case 
of c? > 1, as soon as one can assume that the optimal 
fluctuation of a random potential remains almost uniform 
along t. 
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